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Abstract. This work presents a novel formulation of the ice 
nucleation spectrum, i.e. the function relating the ice crys- 
tal concentration to cloud formation conditions and aerosol 
properties. The new formulation is physically-based and ex- 
plicitly accounts for the dependency of the ice crystal con- 
centration on temperature, supersaturation, cooling rate, and 
particle size, surface area and composition. This is achieved 
by introducing the concepts of ice nucleation coefficient (the 
number of ice germs present in a particle) and nucleation 
probability dispersion function (the distribution of ice nucle- 
ation coefficients within the aerosol population). The new 
formulation is used to generate ice nucleation parameteriza- 
tions for the homogeneous freezing of cloud droplets and the 
heterogeneous deposition ice nucleation on dust and soot ice 
nuclei. For homogeneous freezing, it was found that by in- 
creasing the dispersion in the droplet volume distribution the 
fraction of supercooled droplets in the population increases. 
For heterogeneous ice nucleation the new formulation con- 
sistently describes singular and stochastic behavior within 
a single framework. Using a fundamentally stochastic ap- 
proach, both cooling rate independence and constancy of 
the ice nucleation fraction over time, features typically as- 
sociated with singular behavior, were reproduced. Analysis 
of the temporal dependency of the ice nucleation spectrum 
suggested that experimental methods that measure the ice 
nucleation fraction over few seconds would tend to under- 
estimate the ice nuclei concentration. It is shown that in- 
ferring the aerosol heterogeneous ice nucleation properties 
from measurements of the onset supersaturation and temper- 
ature may carry significant error as the variability in ice nu- 
cleation properties within the aerosol population is not ac- 
counted for. This work provides a simple and rigorous ice 
nucleation framework where theoretical predictions, labora- 


tory measurements and field campaign data can be recon- 
ciled, and that is suitable for application in atmospheric mod- 
eling studies. 


1 Introduction 

Aerosol emissions impact the formation of ice and mixed- 
phase clouds by modifying the background concentration of 
ice-forming particles (e.g. Rosenfeld and Woodley, 2000; 
DeMott et al., 2003a; Pratt et ah, 2009; Prenni et ah, 2009; 
Lohmann and Feichter, 2005; IPCC, 2007; Barahona et ah, 
2010; Hoose et ah, 2010). Satellite retrievals suggest that 
cloud condensation nuclei (CCN) emissions may decrease 
the average size of cloud droplets and “delay” homogeneous 
ice nucleation in convective clouds (Rosenfeld and Wood- 
ley, 2000; Ramanathan et ah, 2001). Ice nuclei (IN) emis- 
sions may lead to enhanced competition between homoge- 
neous and heterogeneous ice nucleation reducing the ice 
crystal concentration in cirrus clouds (Barahona et ah, 2010; 
Karcher et ah, 2006; DeMott et ah, 1994). IN emissions 
can also increase the cloud glaciation temperature in mixed- 
phase clouds leading to enhanced precipitation (Lohmann, 
2002; Lohmann and Diehl, 2006; Diehl et ah, 2007). Ac- 
counting for these effects in atmospheric modeling studies 
requires a link between the ice crystal number concentra- 
tion, cloud formation conditions (i.e. saturation ratio with 
respect to ice. Si, and temperature, T) and the aerosol physic- 
ochemical properties. Such relation is usually termed the 
“ice nucleation spectrum”, N C (S{, T, ni...n) (where p i... n rep- 
resent the moments of the distribution of aerosol number 
concentration) (Barahona and Nenes, 2009; Pruppacher and 
Klett, 1997). Laboratory and field campaign data have been 
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used extensively to generate expressions for N c (Si , T, n | n ) 
(e.g. DeMott et al., 2010; Meyers et al., 1992; Phillips 
et al., 2008; Murray et al., 2011; Welti et al., 2009; Mohler 
et al., 2006; Vali, 1994), however theoretical prediction of 
N C (S{, T, m... n ) remains a challenge. 

Classical nucleation theory (CNT) is often used to calcu- 
late the nucleation rate coefficient, which when integrated 
over the appropriate time scale gives N c (Si, T, /xi ... n ) (e.g. 
Khvorostyanov and Curry, 2004; Liu and Penner, 2005; 
Hoose et al., 2010; Chen et al., 2008; Pruppacher and Klett, 
1997; Barahona and Nenes, 2009). CNT-based models are 
usually associated with the so-called “stochastic hypothe- 
sis” (Pruppacher and Klett, 1997), where all particles in the 
aerosol population are assumed to have the same ice nucle- 
ation probability. Application of this approach to the ho- 
mogeneous freezing of liquid droplets has shown agreement 
with experimental results (e.g. Khvorostyanov and Sassen, 
1998; Liu and Penner, 2005; Koop et al., 2000). However, 
for heterogeneous ice nucleation it provides only a rough 
approximation to the ice nucleation properties of ambient 
aerosol (Marcolli et al., 2007; Murray et al., 2011; Liiond 
et al., 2010) and often requires empirical constraints to re- 
produce observations (Eidhammer et al., 2009; Phillips et al., 
2008; DeMott et al., 2010). 

The assumptions behind CNT however do not require that 
all particles in an aerosol population should have the same 
heterogeneous ice nucleation probability, just that all sur- 
faces of the same composition and structure do (Kashchiev, 
2000). In fact, due to the heterogeneity of ambient parti- 
cles, some variation in the particle surface properties within 
an aerosol population may be expected. This has been rec- 
ognized in recent studies. For example, Marcolli et al. 
(2007) found good agreement between theoretical predic- 
tions and laboratory results when the heterogeneous nucle- 
ation rate coefficient, Thet, was averaged over a distribution 
of contact angles. Working along this line, Niedermeier 
et al. (2011) and Broadley et al. (2012) developed models 
where besides a contact angle distribution, the dependency 
of N c (Si, T,fi i n ) on active site area and the external mix- 
ing of IN were considered. In these models the surface of 
each particle is assumed to be divided into smooth “patches” 
where J\ Kl is locally defined. Using this approach Broadley 
et al. (2012) reproduced experimentally observed dependen- 
cies of N c (Si, T, hi...b) on cooling rate, time and tempera- 
ture in immersion freezing. These models however require 
the knowledge of the area of each surface patch and of the 
intra- and inter-particle distributions of surface composition, 
which may be difficult to constraint, limiting their applica- 
bility in cloud studies. 

A different approach based on the “singular hypothe- 
sis” (Vali, 2008, 1994; Pruppacher and Klett, 1997, and 
references therein), relies on the existence of active sites 
on the surface of the IN where ice nucleation occurs at 
some characteristic T and 6), leading to a time-independent 
N c (Si, T. /i i „ ) (Pruppacher and Klett, 1997). The lack of 


temporal dependency of N c (Si,T, fi\ n ) found in several 
studies (Mohler et al., 2006; Broadley et al., 2012; Connolly 
et al., 2009; Vali, 2008) supports this this view. However a 
theory describing the nature of the active sites that allows the 
prediction of their characteristic T and S\ is yet to be formu- 
lated and N c ( .S'i . T, /i \ n ) is generally fitted to limited exper- 
imental results. Furthermore, ice nucleation does not occur 
instantaneously and the singular hypothesis must be under- 
stood as an asymptotic approximation to rapid ice nucleation 
around the characteristic T and S\. 

Vali (2008, 1994) proposed a hybrid singular-stochastic 
approach for immersion freezing where the IN follow ap- 
proximately singular behavior but there is scatter in the ob- 
served freezing point due to fluctuations in the ice-embryo 
size. This model however does not specify how the main 
freezing point and the width of the dispersion around it 
should be determined. Niedermeier et al. (2011) and 
Broadley et al. (2012) showed that even if the nucleation 
process is fundamentally stochastic, N c (Si, T, /i i n ) may ap- 
proximate singular behavior due to surface composition vari- 
ability. 

There is currently no formulation of the ice nucleation 
spectrum capable of explaining all the observed features of 
ice nucleation within a physically-based framework. This 
work addresses this issue by introducing a novel formulation 
of N C (S\, T, /i] n ) which instead of using a mechanistic ap- 
proximation of the surface structure of the IN, relies on a sta- 
tistical view of the ice nucleation process. The new formula- 
tion is physically-based and explicitly considers the effects of 
supersaturation, temperature, cooling rate, and aerosol size 
and surface area, on the ice crystal number concentration. 
This is achieved by introducing the concepts of ice nucle- 
ation coefficient and ice nucleation probability dispersion. 
The new formulation is applied to generate parameterizations 
of N c (Si, T. //i ...n) for the homogeneous freezing of cloud 
droplets and the heterogeneous nucleation of ice on dust and 
soot IN in the deposition mode. 

2 General theory 

In this section the concepts of ice nucleation coefficient and 
ice nucleation probability dispersion are introduced and used 
to develop a general relation for the aerosol ice nucleation 
fraction, /f (i.e. the fraction of particles that nucleate ice). 
Nucleation theory is then used to link /f to measurable quan- 
tities such as T and 5; and to derive the homogeneous and 
heterogeneous ice nucleation spectra. 

We start by introducing the ice nucleation coefficient, <p, 
i.e., the number of ice germs present within the volume, or 
on the surface of an aerosol particle. The probability of ice 
nucleation, Pf, is heuristically related to <p by (Pruppacher 
and Klett, 1997) 

P f = l-e-v (1) 
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Equation (1) can be extended to describe ice nucleation in an 
aerosol population by finding the weighted average value of 

Pf, 

OO 

f f =l-Jn(<p)e-*d»p (2) 

o 

where n icp) is the distribution function of cp, and describes the 
number fraction of particles in the aerosol population with 
nucleation coefficients between ip and <p + di p. Thus, Eq. (2) 
represents the sum of the contribution of each ip class to the 
overall aerosol ice nucleation fraction. 

The dependency of /f on .S', , T and the physicochemical 
properties of the aerosol population is found by relating ip to 
the “characteristic” nucleation coefficient, ip, chosen to be a 
chemically homogeneous particle with properties (composi- 
tion and size) equal to the average properties of the aerosol 
population at .S, and T . 

Making § = |, Eq. (2) can be written as. 


2.1 Definition of (p and n(§) 

Equation (4) depends only on ip and general expressions for 
ft in terms of ip can be derived without necessarily know- 
ing the structure of the each particle in the population. Thus, 
Eq. (4) holds regardless of the ice nucleation mechanism (ho- 
mogeneous or heterogeneous), or, the interpretation of het- 
erogeneous ice nucleation, e.g. singular vs. stochastic. Its ap- 
plication however requires linking ip to S), T and the aerosol 
properties; this is accomplished below. 

2.1.1 Homogeneous nucleation 

The homogeneous freezing of cloud droplets and deliquesced 
aerosol is generally modeled using the homogeneous nucle- 
ation rate coefficient, /horn- Comparison of Eq. (1) against 
Eq. (7.66) of Pruppacher and Klett (1997) suggests that 

t 

(p — Vp J Thom tit (6) 

0 


oo 

f f =\- (3) 

o 

where n(f)d£ = n((p)d(p has been used (i.e., the number con- 
centration of particles in each (p class must be equal in the f 
and ip spaces). The function n(f ) is termed the “Nucleation 
Probability Dispersion Function” (NPDF), and describes the 
deviation of ip from ip in each particle of the aerosol popula- 
tion. 

It can be noticed immediately from Eq. (3) that 


where n p is the particle volume. Equation (6) can be sim- 
plified by expanding In Thom into its Taylor series around the 
current Si and T, 


In ThomCO = In Thom (Si, T) + - 5)] 

do i 

3 In Thom 

+ J [T(t)-T] + ... (7) 

61 

where the derivatives in Eq. (7) are evaluated at the current 
,S'j and T . Taking the exponential of Eq. (7) gives. 


ff= 1 -N(ip) (4) 

where A f(ip) is the Laplace transform of n(f ). M(ip) and 
n(f) are equivalent representations of the NPDF (Medhi, 
1994). Thus /f contains all the information on the under- 
lying distribution of nucleation coefficients (and vice-versa), 
hence of nucleation probability, in the aerosol population. 

Af(ip) for gamma, sectional and lognormal NPDFs is 
shown in Table 1. For the Lognormal distribution Af (ip) is 
approximated using Af (ip) oc f Q l ^n(%) d£ (Rossberg, 2008) 
which is accurate to within 5 % for a v > 3, when com- 
pared against the direct numerical solution of Eq. (3) (not 
shown). Complex materials may exhibit multimodal NPDFs 
for which n(f) must be modified accordingly (Sect. 3.2.4). 
In the general case of a sectional NPDF calculation of J\f(ip) 
involves a discrete approximation to the Laplace transform 
(Table 1) (Shortle et ah, 2003). 

The ice nucleation spectrum is calculated directly from /f, 

N c (Si, T , /Xi... n ) = NtfmSi, T, Mi...n)] (5) 

where N a is the aerosol number concentration. 


r i /o Ti. / 3 In Thom 9 In Thom \ 

Thom — Thom ('Si, T) exp I + yy, (8) 

where fy = Sft) — .S) and er — T (t) — T are perturbations in 
S, and T, respectively. Introducing the last expression into 
Eq. (6) we obtain. 


/ / 
J Jhomdt = VpJ\i 0m (S[, T) j 


9 In /hom ^ In /horn \ , / 

-isr es ' + -^ eT+ -r m 


using Eq. (9) and making r nuc = 
/ 0 exp ^ j) lT -jy°' n CS', + 31n 3 y om 6r + ■ ■ - j dr', Eq. (6) becomes 


ip — UpThomTiuc (10) 

where r nuc is termed the “nucleation timescale”. In general 
r nu c depends on the predominant nucleation mechanism, T, 
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S; and the cooling rate (Barahona and Nenes, 2008). Intro- 
ducing r nuc in the form of Eq. (10) simplifies the mathemati- 
cal treatment of cp, particularly when S, and T change during 
ice nucleation (Sect. 3.1.1). 

Using as characteristic state a homogeneous droplet with 
the composition and size equal to the mean of the droplet 
population (Sect. 2) we obtain, 

^ fpAiomfnuc (11) 

fp Atom blue 

Equation (11) can be simplified by taking into account that 
ip is dominated by Atom and only linearly dependent on r nuc 
(Barahona and Nenes, 2008). Thus, to a good approximation, 

$ = ( 12 ) 

ftp Atom 

Homogeneous nucleation is a stochastic process (Pruppacher 
and Klett, 1997) where Thom is the same in all droplets that 
have the same composition. If the variation in droplet com- 
position across the population is small (as for example in the 
case of cloud droplets), then 

l>n 

I = = (13) 

i> P 


sites using s p as the characteristic length associated with p as . 
If known, the active site surface area (Pruppacher and Klett, 
1997; Zobrist et ah, 2007) is also suitable to define the contri- 
bution to (p from nucleation on active sites. It must be noticed 
that a mechanistic approximation of (p is not required to de- 
fine the NPDF and it is only used to investigate its physical 
significance. 

The NPDF can be obtained from Eq. (4) by choosing as 
characteristic state a particle with area ,v p and composition 
equal to the bulk mean of the aerosol population. It is as- 
sumed that ice nucleation at the characteristic state occurs 
only through surface adsorption, implying that <p is time- 
dependent. Several reasons justify this assumption. Experi- 
mental studies in immersion freezing have reported temporal 
dependency of N C (S{, T, , n ) (e.g., Broadley et ah, 2012). 
Also robust, time-dependent models (e.g., CNT) exist to de- 
scribe ice nucleation, however only empirical correlations 
(which are typically obtained using aerosol samples of un- 
known heterogeneity) are available for p as . Finally, even if 
nucleation on active sites occurs very rapidly, it is still time- 
dependent. Thus, expressing nucleation on active sites in 
terms Aid (Fletcher, 1969; Pruppacher and Klett, 1997; Mar- 
colli et ah, 2007) would simplify the mathematical descrip- 
tion of (p (Eq. 17). Notice that even if ip is time-dependent, 
the NPDF may be such that on average the temporal depen- 
dency of /f vanishes (Sect. 4.3.2). 

Based on the above description, ip can be written as. 


Thus for constant droplet composition, variability in <p origi- 
nates only from dispersion in n p and the NPDF is determined 
by the droplet number volume distribution (Sect. 3.1). 

2.1.2 Heterogeneous nucleation 


r 

Aid df 

0 

Combining Eqs. (14) and (15) we obtain 



(15) 


The probability that ice nucleates heterogeneously on the sur- 
face of an aerosol particle is usually expressed in terms of 
the heterogeneous nucleation rate coefficient. Aid. or alterna- 
tively, in terms of the surface density of active sites, p as . Aiet 
and p as are locally defined and can vary over the surface of 
each particle as well as among particles in the aerosol popula- 
tion (Niedermeier et ah, 2011; Zobrist et ah, 2007; Broadley 
et ah, 2012). Taking into account that ice nucleation can be 
initiated either by the presence of an active site or by stochas- 
tic adsorption of water molecules (e.g., Vali, 1994), a mech- 
anistic approximation of <p can be written in the form. 


<P = s p 



(14) 


where ,v p is the particle surface area. Equation (14) is ob- 
tained by assuming that the particle surface is divided into 
N “smooth patches” of constant composition, with area .yp 
and local nucleation rate equal to Aiet.k (Niedermeier et ah, 
2011). Equation (14) also accounts for the presence of active 


S P (ais + Ek=l Sp Jo -Aiet.kdf ) 

| > _ p (16) 

Sp JO *^het dt f 

where Aid in Eqs. (15) and (16) is evaluated at the mean 
aerosol properties. Equation (16) can be simplified by as- 
suming that each active site is associated with a patch of sur- 
face area with contact angle close to 0° (Fletcher, 1969), 


(Sk=l 2 fo •Aiet.kdf'') 

% « — — = —f L (17) 

■I p /( ) Aiet df 

where L is the total number of surface area patches in each 
particle, including active sites. Further simplification can 
be achieved by recognizing that ^J/ =1 ^Aiet.k is simply the 
weighted average of Aiet.k over the particle surface, i.e. 


l Aiet.p df 
0 



Aiet.k df 


(18) 
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where /het.p is the “particle effective” nucleation rate coeffi- 
cient. Introducing the last expression into Eq. (17) we obtain. 


% 


./o -Aiei.p dr 
■Sp Jo Thet dC 


(19) 


Equation (19) indicates that if Jhet.p is the same in all parti- 
cles of the aerosol population, then dispersion in ip originates 
only from variation in the particle surface area and the NPDF 
is determined by the number area distribution. 


vapor pressures at T, respectively (Murphy and Koop, 2005). 
Taking into account that for the droplet population .S w ~ 1 
and dSw ^ 0, Eq. (23) can be simplified to 

t T 

J J hom dt' = J^dT' (24) 

0 T 0 

where y — — ^ is the cooling rate, assumed constant during 
nucleation. Using Eq. (8), the last expression can be written 
as 


3 Application 

3.1 Homogeneous freezing of cloud droplets 

Aerosol emissions impact the cloud droplet size distribution 
altering the cloud glaciation temperature and the cold gen- 
eration of precipitation (Rosenfeld and Woodley, 2000; Ra- 
manathan et al., 2001). It is therefore important to determine 
how the droplet volume distribution and the cooling rate play 
a role in determining /f in cloud droplets. Austin et al. (1995) 
suggested that the droplet volume distribution can be ade- 
quately represented by, 

n(§) = (1 + V) + gV~ (1+u)? (20) 

ru + v) 

where £ = ^ and v is related to the relative variance of the 

Up 

droplet volume distribution (varr = [7(v+ 1)] -1 Austin et al., 
1995). Using Table 1 with ip — UpThomTnuc and making a — 
P = (y + 1), we obtain 

M + (21) 

V (1 + V) J 

3.1.1 Nucleation time scale 


r T 

f t aJ Aiom f ( 9 In Thom \ , 0 - 

J hom dt = /exp I — — — er + ...ld T (25) 

0 T 0 

Over the small T interval in which ice nucleation occurs (typ- 
ically about 2 K, Fig. 2) dln d ^° m is almost constant (Barahona 
and Nenes, 2008) and the high order derivatives of In Thom 
vanish. Using this, and taking into account that Thom de- 
creases steeply with T , Eq. (25) can be simplified to 

t T —€j 

f T , / '/horn f / 9 In Jfo om / \ A rj-,r 

/ Aom dr a* / exp I — — — e T j dT (26) 

o T 

where er — T (/) — T . Since dT — der, Eq. (26) can be writ- 

ten as 

/ -€t 

f T i ,/ Thom f ( 9 111 Thom , \ | / 

/ Thom dr « / exp I — — — e T J de T (27) 

o o 

Solving the integral in Eq. (27) and using 
exp 31n 3 ^° m er^ <jC 1 (Barahona and Nenes, 2008) 
we obtain. 


During an experiment at constant Si and T, Eq. (6) is readily 
integrated to give 

ip — fpTjiom Ar ex p (22) 

where Ar ex p is the experimental time of observation. During 
cloud formation r nuc (Eq. 10) must be used instead as T, 5) 
vary during nucleation. 

To develop an expression for r nuc it is advantageous to 
write / ( { Thom dr' in terms of dT/dt and 3 S w /9r, being S w 
the saturation ratio with respect to water, i.e., 

t T q S w 

f Thom dr = f J ^dT'+ f ^PdS' w (23) 

0 T d > ~Sf 

where Tq = 273 K, S^; at = Ps,i/Ps,w is the value of S w at S; = 
1, and p sv , and p s -, are the liquid water and ice saturation 


r 

J Thom dr 
0 


‘/horn 


d In Jhom 

si y 


(28) 


Using the chain rule, where 

Thom is given by the parameterization of Koop et al. (2000). 
From Murphy and Koop (2005), | = -61 32.9/ T 2 

K. Introducing this into Eq. (28) it can be readily seen 
that r nuc = (-“^1®)/)-’ = \ Collect- 

ing terms into Eq. (21) we obtain the final expression for 
homogeneous ice nucleation in cloud droplets, 


A f(<p) = 


fn An 


( 6132.9 y dlnJ hom y 
,om V T 2 dS w ) 

(1 + v) 


■ -(1 + V ) 


(29) 
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3.2 Heterogeneous ice nucleation in the 
deposition mode 

Ice nucleation by direct adsorption of water vapor onto the 
aerosol surface is termed deposition ice nucleation. This pro- 
cess is thought to significantly impact the formation of cirrus 
at low T (Barahona et al., 2010; Hoose et al., 2010; Mohler 
et al., 2006). Weakly wettable species like soot may also 
nucleate ice in this mode at higher T than required for cir- 
rus formation, potentially impacting the formation of mixed- 
phase clouds (Seisel et al., 2005). Currently only CNT pro- 
vides sufficient detail to allow the calculation of /het over a 
wide range of conditions and is therefore used here to derive 
<p for deposition ice nucleation. 

3.2.1 Ice nucleation rate coefficient 


Curry, 2004; Pruppacher and Klett, 1997) however current 
parameterizations depend on largely unconstrained parame- 
ters and such effects are not considered. 

Despite the simplicity of Eq. (30) calculation of /h e t is sub- 
ject to uncertainty. The form of Z h e t, ci, s , and //, depends on 
whether the ice germ grows by direct incorporation of water 
molecules from the vapor phase or they first adsorb onto the 
particle surface and then diffuse to the ice germ. The value 
of ci, s is also strongly influenced by the formation of water 
monolayers on the particle surface (Maattanen et al., 2005; 
Seisel et al., 2005). Maattanen et al. (2005) have shown that 
these factors may introduce up to a factor of two uncertainty 
in / f . 

The steady-state surface concentration is determined by 
equilibrium between the incoming and outgoing molecule 
fluxes (Pruppacher and Klett, 1997), 


Within the framework of CNT, J\ Kl is given by (Pruppacher 
and Klett, 1997; Kashchiev, 2000), 

/het = Zhetci, s /?hexp (3°) 

where Z h e t is the Zeldovich factor, ci s is the surface con- 
centration of water molecules, is the rate at which water 
molecules are incorporated into the critical ice cluster, and k 
is the Boltzmann constant. The remaining term in Eq. (30) 
is the energy of formation of the ice germ, given by (Prup- 
pacher and Klett, 1997) 

Agg = /^cr i/v r 2 (31) 

where cri/ v is the surface tension of ice (106 mJm -2 Prup- 
pacher and Klett, 1997), and r g is the ice germ size given by. 


2u w (Ti/v 

it Tin $ 


(32) 


where v w is the volume of a water molecule (Zobrist et al., 
2007). 

The compatibility parameter, /, in Eq. (31) accounts for 
the reduction in the work of ice formation caused by the het- 
erogeneous surface (Kashchiev, 2000). It is given by 


/ 


1 

-(2 + cost?) (1 — cos 0) 


2 


(33) 


where 6 is the local contact angle between the solid surface 
and the ice germ (Pruppacher and Klett, 1997). In writing 
Eq. (33) it is assumed that the substrate has planar geometry. 
Given the complex geometry of atmospheric particles (e.g. 
Dymarska et al., 2006; Zimmermann et al., 2008) it is not 
clear whether assuming a spherical substrate (or any other 
simple geometry) would represent a better approximation to 
/ than Eq. (33). Therefore curvature effects are not con- 
sidered. Misfit strain can reduce /het (Khvorostyanov and 


«c Pv / A g d \ 

v s \J2jtm w kT \ kT ) 


(34) 


where Agd is the desorption energy of water molecules from 
the particle surface, p v the water vapor pressure, u s the 
molecular frequency of vibration (10 13 s _1 ), and m w the 
mass of a water molecule (2.99 x 10 -26 kg) (Pruppacher and 
Klett, 1997). The mass accommodation coefficient, a c , has 
been introduced into Eq. (34) to account for the low stick- 
ing efficiency of the water vapor molecules onto the particle 
surface (it is assumed that a c = 1.0 for the ice surface) (Prup- 
pacher and Klett, 1997; Seisel et al., 2005). Assuming direct 
water vapor deposition onto the ice germ (which typically 
results in a conservative estimate of /het) gives (Maattanen 
et al., 2005), 


/het — 


«C Pyfw 

V7 m w kT v s 




exp 



(35) 


Equation (35) resembles Eq. (1 1) of Chen et al. (2008), how- 
ever in Eq. (35) / is raised to the —1/2 power instead of 
1 /2; the latter results from neglecting the effect of the solid 
surface on r g (Vehkamaki et al., 2007). The expression of 
Chen et al. may however lead to unphysical behavior as it 
implies that /het -* 0 when f 0 (i.e. ice nucleation would 
be more difficult on easily wettable materials). The authors 
also suggested that Agj is similar to the activation energy of 
water molecules in immersion freezing. However whereas 
increasing Agd leads to low desorption rates increasing ci, s 
(Eq. 34) hence /het. increasing the activation energy leads to 
a decrease in /het (Pruppacher and Klett, 1997). 

Equation (35) can be written in the form 


/ het = A exp 



(36) 


where A = 77 (#)■ T yp ical values ° f ^ 

range between 10 23 and 10 26 m -2 s -1 . 
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3.2.2 Nucleation time scale 

For ice nucleation in the deposition mode the relative varia- 
tion in ,S' W is significantly larger than the relative variation in 
T (Barahona and Nenes, 2009), i.e. y- yy TW' Equa- 
tion (23) can then be approximated as. 


/ *»*'= [ l£< 

0 V 


(37) 


Using yf- a s M (Barahona and Nenes, 2009), where u is 
the vertical velocity and a s — — yy, g is the accel- 

eration of gravity, A// s the latent heat of sublimation of ice 
(Pruppacher and Klett, 1997), and c p and m a are the heat 
capacity and molecular mass of air, respectively, Eq. (37) be- 
comes. 


l o w 

/ /hetdf — / ./her d 

J a s u J 

0 Sisat 


(38) 


Replacing 7 hom for /het into Eq. (8) and introducing the result 
into Eq. (38) we obtain 


0 y 

To calculate is convenient to write r g in terms of S v , 

Making S w — Si (p s ,i/.Ps,w). replacing into Eq. ( 32 ) and re 
arranging we obtain. 


1 + 


s„-i 


(43) 


ln(Ps, i/Ps,w) 


Where r°=r„ L and ln(S w ) S w — 1 has been used. 
Taking the derivative with respect to S w in Eqs. (31) and (43), 
replacing into Eq. (42), and rearranging we obtain. 


3 In /he t 

3 S w 


/< 


(l+ Sw-l V 

V ln(/> s ,i/Ps,w)/ 


(44) 


„o _ 4 n {‘if 


where n° = 7 is the number of water molecules in the 

ice germ at water saturation. The second term in the denom- 
inator of Eq. (44) is typically negligible compared to 1, and 
Eq. (44) can be approximated as 


3 In /het 

9S W 


fnl 


(45) 


Introducing the last result into Eq. (41) gives, 


[ /het dr' = — f exp 

J a s u J 


/ 3 In / h et 

V dS* 




~ e S w 


d S' 


(39) 


where e^ w = S w (t) — S w . For small variation in 5 W , 81 3 n ^ hel is 
almost constant (Barahona and Nenes, 2009) and the high or- 
der derivatives of In /h et vanish. With this Eq. (39) becomes. 


t o 

f Jh « dt ' =: tu / exp (^ W )< < 46 > 

0 -es w 

Integrating the last expression, and since typically 
exp /n°es w ^ 1, Eq. (15) can be written as 


f /het dr' = — f exp f 

J a s u J \ 


( 9 In / het 


9S W 




d S' 


(40) 


Sj* at 


Taking into account that /het is a steep function of S w 
(Sect. 4.3) and that /h e t = 0 for S w — S^ at the lower limit 
of integration in Eq. (40) can be replaced by S w — €s w ■ Using 
this, and making dS’w = des w into Eq. (40) we obtain 


/ Jh » d, ' = ^ / exp (^ ei - ,d6 /- 


( 41 ) 


~ e S„ 


Further simplification can be achieved by noticing that most 
of the variability in /h e t comes from the exponential term in 
Eq. (36). Therefore to a good approximation we can write 


(p — -Vp /het I 


1 


o J — s p /het Flue 


(47) 


^a s ufn^ 
where r nuc = (a & ufn° g )~ x . 

3.2.3 The nature of /t(§) in deposition ice nucleation 

The mechanistic approximation of <p presented in Sect. 2. 1 .2 
can be used to investigate the functional form of the NPDF. 
Using Eq. (47) into Eq. (19), and since most of the variability 
in (p comes from /het.p- we can write, 

| (48) 

■Sp/het 

Using Eqs. (31), (35) and (48), and taking into account that 
most of the variation in /het.p results from variation in the 
exponential term of Eq. (36), we obtain for the NPDF 


9 In /het 1 9Ag g 3 r g 

35’w kT 3 r g dS w 


(42) 


(Z) = „ HZe-AeiWpy-m) 


(49) 
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Table 1. A f (<p) and n(f) for gamma, lognormal and sectional NPDFs; <p is the characteristic ice nucleation coefficient defined in Sect. 2.1, 
T is the gamma function, and / a £ represents the fraction of the aerosol population with cumulative ice nucleation probability below 1 — 

exp(— %pip). The lognormal Af(ip) is approximated using Af (<p) oc d£ (Rossberg, 2008). 


Distribution 

n(M) 

A/( <p) 

Gamma 

r(ce) 5 e 

( _j_\ a 

Lognormal 

In 2 (if) 

2<tJ 
e v 

1 erfc ( Xn{ ^ \ 



Sectional 

fa,k~ Ja,k— 1 

Y^Pi e ~^‘ where Pi = 1 1,1 (§) d ? 

i 


where 9 p is the “particle effective” contact angle associated 
with /het.p (Eq. 19), 9 is the contact angle associated with 
ip, and Ag° = Agg| 0=lgoo . Expanding f{9 p ) into its Taylor 
series around 9 and neglecting the high order derivatives of 
/ gives. 




(50) 


which can be written as 


n(g) = n( S ±e« e r-V 


(51) 


where c = -Ag° % 


If all particles in the population have 


the same surface area, the NPDF reduces to 


n(f) = n {e c(e p- 0) ^ (52) 

Equation (52) suggests that «(§) for heterogeneous ice 
nucleation follows a lognormal distribution. If Op — 0 fol- 
lows a normal distribution with variance og (i.e. the parti- 
cle surface composition is randomly determined), then £ = 
exp (c(9 p —9)) follows a lognormal distribution with variance 
c 2 a £ (this can be shown by making lr '^ = ( 6 p — 9)). At the 
limit of constant 9 p the number area distribution (which is 
usually lognormally distributed, Seinfeld and Pandis, 1998) 
determines n(§) (Eq. 19). These two limits of variability in 
<p suggest a lognormal NPDF. Equation (52) also indicates 
that 9 is approximately equal to the mean 9 p of the aerosol 
population. 

3.2.4 Final form for deposition ice nucleation 

Combining the expressions of Table 1 with Eqs. (4) and (47), 
a concise expression for M {ip) can be written in the form. 


J\f (cp) — -erfc 


In ( 

V°w'ig/ 

^~7T G(p 


(53) 


where / and 7het are calculated at the mean aerosol proper- 
ties, and a v is termed the “ice nucleation dispersion coeffi- 
cient”. Equation (53) shows that given 5;, T and u, ff de- 
pends only on s p , a ([! , and the characteristic contact angle, 9 
(associated with ip). Complex aerosol mixtures may exhibit 
multimodal NPDFs. Equation (53) can be readily extended 
to such cases by making. 


j M 

N(<p) — ~ w ' erfc 


i=i 


\“s ufin°J 
y/lt &(p.i 


(54) 


where Wi is the weight of the i-th lognormal mode in n(f), 
and M is the total number of modes. 


4 Results and discussion 
4.1 General behavior of A f(<p) 

The probability of observing a nucleation event in an aerosol 
population increases with ip, i.e. ff is always a monotonic 
function of ip. Figure 1 shows /f for lognormal, gamma, 
and sectional NPDFs. The width of each NPDF (e.g. a, 
and the standard deviation, er, respectively) has been varied 
over a wide interval and the distribution mean set to f = 1 
(Eq. 3). To illustrate the discrete transform approximation 
method (Table 1, Shortle et ah, 2003), the sectional distribu- 
tion in Fig. 1 (top panels) has been derived from a normal 
distribution using 20 bins. This method produces a contin- 
uous, smooth increasing ff, however using only 20 bins re- 
sults in smaller sensitivity to a than the analytical transform 
of n(f ) (not shown). 

The form and width of the NPDF strongly influence ff. 
As a consequence of the support of the normal distribution 
between — oo and oo, nucleation spectra derived from the 
normal NPDF tend to shift right as o decreases. Indeed, if <p 
is normally distributed then a finite probability of cp < 0 ex- 
ist, leading to unphysical behavior. As the gamma and log- 
normal NPDFs have support only on the positive real axis. 
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Fig. 1. Ice nucleation probability dispersion functions (left panels) and corresponding ice nucleation fraction (right panels). Shown are 
sectional (based on a normal distribution using 20 bins, top panels), gamma (middle panels) and lognormal (bottom panels) distributions. 





Fig. 2. Homogeneous ice nucleation spectra for the freezing of cloud droplets. Unless otherwise specified y = — ^ = 1 Kmin 1 , Dp = 20 
pm, and varr = 0.25. 
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Fig. 3. Heterogeneous ice nucleation spectra in the deposition 
[5°, 10°, 15°, 23°] (from left to right). Unless otherwise specified u = 

values of cp < 0 are not accessible to N (ip). Thus, the inflec- 
tion point in the /f curve (Fig. 1) for these NPDFs is always 
around ip = 1 . This means that ice crystal formation becomes 
significant around the conditions for which there is, on aver- 
age, one ice germ per particle in the population. <p = 1 also 
defines characteristic values for Thet (or Thom) and p as . A 
similar conclusion was reached by Khvorostyanov and Curry 
(2004) using the empirical constraint of s p /het % Is -1 . The 
theory presented here generalizes this picture and suggests 
that the characteristic 7het (at ip = 1) is found at /f ~ 0.5 in- 
stead of the typical /f ~ 0.01 used in experimental studies. 

The lognormal and gamma distributions display distinc- 
tive behavior to the variation in the width of the NPDF. For 
the gamma distribution the slope of /f remains almost con- 
stant for a > 0.5 but decreases steeply for a < 0.5, which 
is explained by the larger probability of finding small (p at 
lower a (Fig. 1, left panels). For the lognormal distribution a 
larger leads to a larger probability of both, small and large 
<p. The former limits the value of fi as — > oc whereas the 
latter increases /f at low ip. These features may explain the 
low /f typically observed in ambient aerosol (e.g. Eidham- 
mer et al., 2009; Phillips et al., 2008) and the diversity of ice 
nucleation thresholds observed in experimental studies. This 
is further analyzed in Sect. 4.3.3. 



S, 

mode. Lines are grouped by characteristic contact angle, 9 — 

0.1 ms -1 , Jp = 5 pm 3 , a v = 5, and T = 225 K. 

4.2 Homogeneous nucleation 

Figure 2 shows the sensitivity of the nucleation spectrum for 
homogeneous freezing to D p , y and varr (Eq. 21). Although 
Eq. (29) is written in terms of E p , the equivalent droplet size, 
Dp, is most often used in experimental studies; they are sim- 
ply related by D p = ■ As ip scales with y and D p 

their variation tends to “shift” /f along the T axis. Decreas- 
ing D p from 40 pm to 10 pm decreases 7 onsct (defined at 
f{ ~ 0.01) from 236 K to 233 K. Variation in y has the op- 
posite effect although with a weaker impact: a factor of 10 3 
increase in y only results in about 3 K lower 7 onsct (explained 
by the decrease in r nuc as y increases, Eq. 27). Thus cooling 
rate variation may have a limited role in ice crystal produc- 
tion by homogeneous freezing of cloud droplets (although it 
has an important effect in the homogeneous freezing of liquid 
aerosol particles (Barahona and Nenes, 2008)). 

Homogeneous freezing may be strongly impacted by dis- 
persion in the droplet size distribution. For varr < 0.4 the 
slope of /f remains almost constant and nearly all droplets 
freeze within 2-3 K of 7 onsct . However as varr increases 
beyond 0.4, fi “stretches” over a wider T interval so that 
at varr = 0.7 some droplets freeze at T as low as 228 K. 
This is explained by the decrease in a as varr increases (i.e., 
varroc T) resulting in a larger fraction of small droplets in 
the population (Sect. 4.1). In other words since T onset does 
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Fig. 4. Ice nucleation fraction vs. observation time at T = 220 K, for IN with 8 = 6° (left panels) and 9 = 23° (right panels), (a) Constant 
5; = 1.023. (b) Constant 5; = 1.345. (c-d) 5; is increased over 1 s up to the point where /f = 0.001 and then maintained constant afterwards, 
(e — F) As (c) and (d) but 5) is increased by 0.001 after Af ex p = 100 s. The particle mean surface area was set to Jp = 10 pm 2 . 


not change with increasing varr the stretching effect becomes 
only obvious if the whole /f curve is considered. This pre- 
viously unidentified behavior of homogeneous freezing may 
have important implications for the development of convec- 
tive clouds as aerosol emissions and entrainment tend to 
broaden the droplet size distribution (Pruppacher and Klett, 
1997). 

4.3 Heterogeneous ice nucleation in the 
deposition mode 

4.3.1 General features 

Figure 3 shows representative profiles of /f for heteroge- 
neous ice nucleation in the deposition mode. In general /f 
is primarily a function of .S', and 9, and to a lower extent, 
of T, .Vp, u and a v . For 9 < 10°, /f tends to be very steep 


and almost insensitive to T, ,Vp, u and a v . Therefore a sin- 
gle, constant Sunset can be assigned to aerosol populations 
with low 9. This indicates that highly efficient IN tend to 
display features typically associated with singular behavior 
(even though a time-dependent approach is used to calculate 
ip). However this may not imply a lack of temporal depen- 
dency of /f at constant .S’, and T (Broadley et al., 2012) and 
is further analyzed in Sect. 4.3.2. 

IN with high 9 are strongly influenced by T, s p and u. 
In this regime /f is also impacted by so that variation in 
the width of the NPDF caused by repeated freezing cycles 
and by dispersion in surface composition may affect S) onset- 
These characteristics imply that weak IN display features 
commonly associated with the stochastic hypothesis. Thus 
within the context of the theory presented here, stochastic 
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c 

*“ 7 , onset 


Fig. 5. Correlation between the mean contact angle and 5; at the 
nucleation onset (ff = 0.01). T was varied between 200 and 250 K, 
u between 0.01 and 1 ms -1 , and ip between 0.1 and 100pm 2 . 


and singular features in the ice nucleation spectrum can be 
understood as limits of variability in the NPDF. 

For the highest 9 used in Fig. 3 (23°) decreasing T from 
250K to 210K leads to an increase of 0.15 in Si, 0 nset (de- 
fined at f{ Ri 0.01). This is explained by the increase in r g 
and Ag g as T decreases lowering /het (Eq. 30). Increasing 
u from 10 -3 to 10ms _1 increases Sunset by about 0.1 at 
9 = 23° (Fig. 3, bottom left panels). This is explained by 
the decrease in r nuc as u increases (Eq. 47) so that higher 
/het is required to reach (p = 1 . Although the sensitivity of 
Si, onset to u is significant, achieving a large variation in u in 
experimental studies may be technically challenging and the 
effect of variation in the cooling rate on ff in deposition ice 
nucleation may be difficult to observe. 

Si, onset decreases by about 0.08 when .? p increases from 
0. 1 pm 2 to 100 pm 2 . Although this is consistent with existing 
studies (Kanji et al., 2008; Welti et al., 2009) direct compar- 
ison against experimental results is difficult as the number 
area distribution is usually not reported. The range of ,v p con- 
sidered in Fig. 3 is however typically observed in ambient 
aerosol (Eastwood et al., 2008; Dymarska et al., 2006; Kanji 
et al., 2008) and the effect of s p on ff may be readily observ- 
able in the laboratory. 

Increasing a v tends to decrease S M)nsct as larger a v leads 
to larger probability of finding high ip (Sect. 4.1). This effect 
is evident for 9 > 10° and indicates that the relation between 
Si, onset and 6 (Fig. 5) can be strongly impacted by a^. When 
<j<p —> 0 (i.e. a very narrow NPDF, Eq. 1) 0 and Sj, onset are re- 
lated by an exponential function weakly dependent on 7', ij, 
and u, and close to the solutions of /het = - l!l (Q,99) . A sim- 

Sp^nuc 

ilar functional dependency was found by Wang and Knopf 
(2011) from measurements of S;. onset, which was interpreted 


as a physical dependency of 6 on ,S). Figure 5 instead indi- 
cates that rather than physical dependency the empirical fit of 
Wang and Knopf represents correlation. Indeed, for cr (p > l, 
different combinations of T , ,s p and u result in a family of 
relations between 9 and S;, onset- 

Figure 5 suggests that the error in 9 from neglecting vari- 
ability in surface properties increases with 5;, onset and 
Determining Sj, me d (he., S\ at ff — 0.5) instead of Sunset 
may reduce this error (Fig. 3), however only measurement 
of the ff vs. S) curve provides information on rr^, and 
would indicate the existence of multiple modes in the NPDF 
(Sect. 4.3.3). 


4.3.2 Temporal effects on deposition ice nucleation 


Figure 4 shows ff for efficient and weak IN (0 = 6° and 
0 = 23°, respectively) as a function of Af exp for different 
values of S; and constant T. For these conditions Eq. (53) 
is simplified as N(ip) = i>crfc ln( ’ ;p ^ 1 j Xr ‘ : ’ Lp> Maintaining 
the efficient IN at S, = 1.023 (Fig. 4a) results in significant 
ff after a few seconds. For low surface variability (a v < 5) 
ff keeps increasing overtime whereas foro^ > 10 it remains 
constant after Af ex p = 10 s. This behavior is explained by the 
larger increase in < p required to increase f f as o v increases 
(Eqs. 37 and 53). It is also consistent with the work of Nie- 
dermeier et al. (2011) who concluded that increasing vari- 
ability in surface properties weakens temporal effects. The 
temporal behavior of ff is strongly influenced by the value 
of Si set during an experiment. Thus maintaining the IN 
population with 9 — 23° at Si = 1.023 results in ff — 0 over 
the whole time interval considered (not shown), and only for 
Si — 1.345 (Fig. 4b) is the temporal dependency of ff evi- 
dent. 

The value Si at which temporal effects on ff become evi- 
dent may not be known in advance. Thus a more likely Si tra- 
jectory in experimental studies may increase Si up to a point 
where ff = 0.001, and then maintain it constant afterwards 
(Figs. 4c and 4d) (e.g. Broadley et al., 2012). For this case 
IN populations with low result in high ff over the time of 
the experiment since they start nucleation at higher 5) than 
IN with higher This behavior is opposite to the cases at 
constant Si (Figs. 4a and 4b) where low resulted in low ff, 
and indicates that the 5) trajectory followed in an experiment 
strongly impacts the observed ff. As before, high cr v , leads to 
weak temporal dependency of ff. 

High variability in surface properties may not always re- 
sult in a weak temporal dependency of ff. In Fig. 4e and 
4f the same conditions as in Fig. 4c and 4d are maintained, 
but 5) is increased by 0.1% after Af cxp = 100 s. For 9=6° 
(Fig. 4e) this perturbation in S; rapidly increases ff which 
remains constant afterwards. For 9 = 23° (Figs. 4f) a simi- 
lar increase in ff occurs after the increase in 5), however for 
o,,. > 1 0 temporal effects become more evident after the per- 
turbation. The striking feature of Fig. 4f is that it is the IN 
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Fig. 6. Heterogeneous ice nucleation spectra in the deposition mode for Kaolinite (KAO) and Montmorillonite (MONT) (Table 2). 


with high Gy that display significant temporal dependency, 
contrary to the behavior depicted in Figs. 4a-d. Thus the 
role of surface variability on the temporal dependency of /f 
depends strongly on the Si trajectory followed. This also in- 
dicates that for IN with multimodal NPDF (Sect. 4.3.3) dif- 
ferent fractions of the aerosol may display different temporal 
behavior. 

Figure 4 suggests that experimental apparatus with A f exp 
of a few seconds would tend to underestimate /f because 
only the fraction of the IN that nucleates ice rapidly at the 
experimentally set Si would be accounted for. Whether this 
bias would affect IN with high or low depends on the Si 
trajectory followed. The magnitude of the underestimation 
depends on a ^ but in some cases it may be of several or- 
ders of magnitude. This indicates that measurements of IN 
in cloud chamber experiments with typical residence times 
of about 10-20 s, would tend to underestimate the IN con- 
centration if the temporal dependency of N c (Si, T, n \ ... n ) is 
not accounted for. 

Comparison between Figs. 3 and 4 suggests that high sen- 
sitivity of /f to u does not necessarily imply significant tem- 
poral dependency of /f. This is because the effect of u on 
/f is controlled by ^ whereas the temporal dependency of 
/f depends on and s p /het- Thus /f for efficient IN may 
be insensitive to variation in u but still exhibit temporal de- 


pendency. An vice versa, the /f of weak IN may not exhibit 
temporal dependency even if it is susceptible to variation in 
u (Fig. 3). This is consistent with the study of Broadley et al. 
(2012) who found that cooling rate and temporal effects on ft 
are not correlated. Moreover, since the lack of dependency of 
/f on u and time can be reproduced using a time-dependent, 
stochastic approach, they cannot be considered prove of the 
singular hypothesis. 

4.3.3 Parameterization of the ice nucleation spectra of 
dust and soot 

Among the different IN species present in the atmosphere, 
soot and dust have been identified to play a significant role in 
the formation of cirrus and mixed-phase clouds (e.g. DeMott 
et al., 2003a, b; Gayet et al., 2004). The theory presented 
in Sects. 2 and 3 can be used to analyze empirical ice nu- 
cleation spectra and generate parameterizations for usage in 
atmospheric models, as follows. 

Dust 

Several authors have studied the ice nucleation properties of 
dust (e.g. Eastwood et al., 2008; Zimmermann et al., 2008, 
and references therein). Although most studies focus on 
measuring .S’ L onset, the works of Mohler et al. (2006, M06) 
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Table 2. Parameters of the nucleation spectra derived from the data reported by Mohler et al. (2006, M06), Welti et al. (2009. W09). and 
Crawford et al. (2011, Cl 1); is assumed to be the same in all modes. 


Species 

Reference 

ip (pm 2 ) 

T (K) 

Mode 1 
U’l h (°) 

Mode 2 
«’2 02(°) 

Mode 3 

W 3 e 3 (°) 

&(p 

ATD 

M05(IN02-148) 

0.64 

223.3 

0.22 

8.7 

0.34 

11.5 

0.44 

16.9 

12.4 

ATD 

M06(IN02-148) 

0.64 

223.4 

0.21 

10.9 

0.29 

14.4 

0.50 

18.0 

29.8 

ATD 

M06(IN02-149) 

0.64 

223.5 

0.17 

6.7 

0.29 

10.5 

0.54 

14.2 

15.8 

ATD 

M06(IN03-05) 

0.64 

210.8 

0.21 

10.5 

0.30 

12.6 

0.49 

14.7 

8.1 

ATD 

M06(IN03-06) 

0.64 

210.6 

0 


0.41 

11.7 

0.59 

13.1 

5.9 

ATD 

M06(IN03-067) 

0.64 

210.5 

0.02 

10.2 

0.39 

11.6 

0.59 

13.7 

9.3 

SD 

M06(IN04-40) 

0.84 

224 

0 


0.16 

22.9 

0.84 

33.6 

23.5 

SD 

M06(IN04-44) 

0.84 

213.1 

0.27 

9.8 

0.05 

14.0 

0.68 

19.6 

30.0 

SD 

M06(IN03-11) 

0.84 

210 

0.21 

11.2 

0.27 

14.3 

0.49 

17.7 

25.9 

MONT 

W09(200 nm) 

0.25 

233 

0 


0.08 

16.1 

0.92 

21.1 

30.0 

MONT 

W09(200 nm) 

0.25 

228 

0 


0.06 

13.3 

0.94 

16.4 

15.9 

MONT 

W09(200 nm) 

0.25 

223 

0 


0.10 

16.0 

0.90 

19.8 

20.8 

MONT 

W09(200 nm) 

0.25 

218 

0 


0.11 

18.5 

0.89 

23.7 

22.0 

MONT 

W09(400 nm) 

1.00 

233 

0 


0.05 

13.4 

0.95 

16.6 

25.6 

MONT 

W09(400 nm) 

1.00 

228 

0 


0.10 

14.2 

0.90 

18.5 

20.0 

MONT 

W09(400 nm) 

1.00 

223 

0 


0.08 

15.8 

0.92 

20.1 

19.9 

MONT 

W09(400 nm) 

1.00 

218 

0 


0.07 

15.2 

0.92 

21.7 

25.0 

MONT 

W09(800 nm) 

4.01 

233 

0 


0.13 

12.1 

0.87 

14.5 

16.5 

MONT 

W09(800 nm) 

4.01 

228 

0 


0.12 

12.8 

0.88 

15.2 

14.3 

MONT 

W09(800 nm) 

4.01 

223 

0 


0.14 

15.0 

0.86 

18.6 

16.4 

MONT 

W09(800 nm) 

4.01 

218 

0 


0.04 

17.8 

0.96 

22.6 

18.3 

KAO 

W09(200 nm) 

0.25 

233 

0 


0.08 

18.4 

0.92 

22.9 

21.8 

KAO 

W09(200 nm) 

0.25 

228 

0 


0.12 

16.1 

0.88 

21.9 

21.0 

KAO 

W09(200 nm) 

0.25 

223 

0 


0.15 

15.8 

0.85 

21.8 

19.5 

KAO 

W09(200 nm) 

0.25 

218 

0 


0.14 

17.2 

0.86 

22.5 

19.6 

KAO 

W09(400 nm) 

1.00 

233 

0 


0.06 

15.1 

0.94 

20.9 

29.8 

KAO 

W09(400 nm) 

1.00 

228 

0 


0.16 

14.8 

0.84 

19.5 

30.0 

KAO 

W09(400 nm) 

1.00 

223 

0 


0.13 

16.0 

0.87 

20.6 

30.0 

KAO 

W09(400 nm) 

1.00 

218 

0 


0.08 

17.4 

0.92 

21.2 

19.2 

KAO 

W09(800 nm) 

4.01 

233 

0 


0.15 

13.6 

0.85 

19.1 

30.0 

KAO 

W09(800 nm) 

4.01 

228 

0 


0.15 

13.4 

0.85 

17.2 

21.2 

KAO 

W09(800 nm) 

4.01 

223 

0 


0.10 

14.8 

0.90 

18.3 

17.6 

KAO 

W09(800 nm) 

4.01 

218 

0 


0.1 

17.8 

0.93 

21.8 

23.2 

OC30 

C1KIN11-20) 

159.4 

225.1 

1 

32.3 





29.7 

OC30 

C1KIN11-21) 

159.4 

224.8 

1 

35.5 





32.6 

OC30 

C1KIN11-22) 

159.4 

224.5 

1 

34.6 





22.7 


Table 3. Summary of average parameters (standard deviation) for the nucleation spectra of Table 2. 


Species 

Mode 1 

m hn 

Mode 2 

U>2 #2(°) 

w 3 

Mode 3 

h (°) 

G(p 

ATD 

0.17(0.08) 

9.4 (1.7) 

0.34(0.05) 

12.0(1.3) 

0.50 

15.1(1.9) 

13.6(8.7) 

SD 

0.24 (0.05) 

10.5(1.0) 

0.16(0.11) 

17.1(5.0) 

0.60 

23.6(8.6) 

26.1(5.1) 

MONT 



0.09(0.03) 

15.3 (2.0) 

0.91 

19.1(2.9) 

20.4 (4.6) 

KAO 



0.14(0.07) 

16.0(1.8) 

0.86 

21.4(3.6) 

25.5 (4.7) 

ALL DUST 

0.19(0.08) 

9.7 (1.6) 

0.15 (0.10) 

15.1(2.6) 

0.66 

19.3(3.9) 

20.8 (7.8) 

OC30 

1 

34.2 (1.7) 





28.3 (5.0) 
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Table 4. Suggested parameters for the parameterization of depo- 
sition ice nucleation on “generic dust” and soot with low organic 
carbon content (OC30). 


Species 

Mode 1 

Mode 2 

Mode 3 

&(p 


«’l 

C(°) 

w 2 

e 2 (°) 

w 3 

%(°) 


DUST 

0.20 

9.9 

0.25 

15.1 

0.55 

19.8 

21.4 

SOOT 

1 

34.2 





28.3 


and Welti et al. (2009, W09) have reported empirical ex- 
pressions for the ice nucleation spectra of dust in the depo- 
sition mode. M06 studied the ice nucleation properties of 
Arizona Test Dust (ATD), Saharan Dust (SD) and Arizona 
Dust (AD) using the AIDA cloud chamber (Mohler et al., 
2003) and reported exponential fits to f{. W09 studied the ice 
nucleation properties of several size-selected dust species us- 
ing the Zurich Ice Nucleation Chamber (ZINC, Stetzer et al., 
2008) and reported sigmoidal fits to /f for Kaolinite (KAO) 
and Montmorillonite (MONT). 

The expressions reported by M06 and W09 are used to de- 
rive the parameters of Af(y>) for the different dust species. 
Basically, 9 (associated with ip) and a v in Eqs. (53) and (54) 
are found for each of the /f curves reported by M06 and W09 
(Table 2). From Seisel et al. (2005), Agd = 6.5 x 10 -20 I and 
a c — 6.3 x 10 -2 . To reduce the parameter space it is assumed 
that a v is the same for all modes. To minimize the effect of 
condensation, only experiments with S w < 0.9 are used. For 

M06 it is assumed that Ar exp = where ^ is the mea- 

sured rate of change of 5; (Mohler et al., 2006). For the W09 
data it is assumed that Af exp = 12 s. M09 reported the value 
of ,?p used in their experiments however W09 reported the 
mean equivalent aerosol size. For the latter, ,y p is estimated 
assuming spherical particles and a lognormal size distribu- 
tion with geometric mean dispersion of 2.5 (cf. Welti et al., 
2009, Fig. 3). The values of 0;, a v and u>, are listed in Table 2 
and the resulting /f vs. Si curves are shown in Figs. 6 and 7. 
Since the number area distribution was not measured in nei- 
ther M09 nor W09, area dispersion and surface composition 
effects are assumed to be combined into rr (/J . 

The NPDF of KAO and MONT is well represented using 
bimodal lognormal distributions (Fig. 6, Table 2). Since ATD 
and SD are mixtures of different materials, three lognormal 
modes are required to represent their NPDFs (Fig. 7, Table 
2). It is however remarkable that such complexity can be 
described using a few parameters. 

Consistent with W09, the spectra derived for the KAO 
sample with D eq = 800 nm tends to reach higher f{ than for 
Z) e q = 200nm (Fig. 6). This is however not the case for 
MONT, which may be a consequence of the slightly lower 
t7(p in MONT than in KAO (Table 3) i.e. high /f is already 
reached in the MONT sample with D eq — 200 nm. However, 
measurement of the number area distribution is required to 
better constraint a v and understand such behavior. 


The average w,, a v and (:), of each mode for the differ- 
ent dust species of Table 2 is presented in Table 3. All dust 
species exhibit modes in the NPDF centered around 9, ~ 15° 
and 9j ~ 20° (although they are centered at slightly lower 
0j in ATD, Table 3). ATD and SD exhibit a third mode at 
9 10°. The standard deviation of 9j among all dust species 

is typically between 10% and 20% of the average 9j. This 
consistency suggests that the NPDF of very different dust 
species may be parameterized using a linear combination of 
a few lognormal modes. In principle w, for each mode would 
be linked to the dust source. Such information is however not 
readily available in most models. A “generic dust” parame- 
terization is therefore formulated by averaging uq and w 2 
between ATD and SD (with w 2 = 1 — ttq — w 2 )- The NPDF 
of dust is assumed to follow a three-modal lognormal dis- 
tribution with 9j and a v averaged among all the species of 
Table 3. The parameterization is summarized in Table 4 and 
shown in Fig. 9. 

Soot 

A parameterization for the ice nucleation spectrum of soot 
can also be derived from experimental data. Few studies 
however focus on ice nucleation on soot particles (Dymarska 
et al., 2006; Mohler et al., 2005; DeMott et al., 1999; Gor- 
bunov et al., 2001). Crawford et al. (2011, Cll) studied 
the nucleation properties of soot derived from the combus- 
tion of propane. It was found that ice nucleated in the de- 
position mode on uncoated soot with organic carbon con- 
tent below 30% (OC30). Crawford et al. (2011) reported 
values of 5; at fi = 0.1% and /f = 1 %, which are used to 
constraint Eq. (53). The particle mean surface area, J p , is 
estimated using bulk surface area and density of 32 m 2 g ~ 1 
and 1600 Kg m , respectively (Popovitcheva et al., 2000), 
and assuming a lognormal size distribution with geometric 
mean diameter of 250 nm and geometric dispersion equal 
to 2 (Crawford et al., 2011). From Seisel et al. (2005), 
Agd = 4.4 x 10 -20 J and a c — 4.7 x 10 -2 . From the Cl 1 
data, Af exp = 80 s, corresponding to a cooling rate of about 
2 Kmin -1 . 

Figure 8 shows the derived nucleation spectra for OC30 
soot. 6 ranges between 32° and 36° and a v between 22 and 
32. The latter is slightly larger than the dust average cs v . 
Due to the fractal characteristics of soot (Dymarska et al., 
2006; Gorbunov et al., 2001), it is possible that surface area 
dispersion plays a more significant role in determining the 
NPDF than for dust. The limiting fi (at 5) = 1 .7) is about 
25 %, in good agreement with the Cll data. The average 9 
and a v among the spectra of Table 2 (OC30 cases) are used 
to parameterize fi (Table 4). The resulting parameterization 
is shown in Fig. 9; as it is based on limited data its usage is 
recommended only for exploratory studies. 
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Fig. 7. Heterogeneous ice nucleation spectra in deposition mode for 
Arizona Test Dust (ATD) and Saharan dust (SD) (Table 2). Lines 
with the same color indicate runs using different aerosol samples of 
the same material (Mohler et al., 2006). 



S, 

Fig. 8. Heterogeneous ice nucleation spectra in the deposition mode 
for OC30 soot (Table 2). 


5 Summary and conclusions 

A novel formulation of the ice nucleation spectrum for ho- 
mogeneous and heterogeneous ice nucleation was developed. 
This was accomplished by introducing the concepts of ice 
nucleation coefficient and ice nucleation probability disper- 
sion. It was shown that the NDPF and the aerosol ice nu- 
cleation spectra are simply related by Laplace transforma- 
tion. The new formulation accounts for the dependency of 
N C (S[, T, [i i „) on particle size, T , S), u, and, in the case of 
heterogeneous ice nucleation, on the distribution of particle 
area and surface composition. It was applied to the homoge- 
neous freezing of cloud droplets and the heterogeneous ice 
nucleation on IN in the deposition mode. For the latter, pa- 
rameterizations of deposition ice nucleation on dust and soot 
for usage in atmospheric models were developed. 



s, 


Fig. 9. Parameterization of the heterogeneous ice nucleation spectra 
for deposition ice nucleation on “generic dust” (DUST) and soot 
with low organic carbon content (SOOT) (Table 4). The shaded 
area corresponds to T = [200-250] K, u = [0.01 — 1] ms -1 , and 
s p = [0.1-100] pm 2 . 


It was shown that the NPDF for the homogeneous freez- 
ing of cloud droplets is determined by the droplet volume 
distribution and well represented by a gamma distribution. 
Analysis of AdS;, T, /i i... n ) for this case showed that varia- 
tion in Dp and y may result in about 2-3 K variation in 7’ onsct . 
Dispersion in the droplet volume distribution however does 
not impact T a nse t but “stretches” /f over a wider T interval so 
that some droplets in the population may freeze at T as low 
as 228 K. This effect is significant for varr > 0.4 (Eq. 20) 
and may have important implications for the cold generation 
of precipitation as CCN emissions and entrainment tend to 
broaden the droplet size distribution. 

The theory presented here suggests that the NPDF for 
heterogeneous ice nucleation on simple materials follows a 
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lognormal distribution whereas the NPDF of complex mate- 
rials can be adequately represented by a sum of lognormal 
modes. It was also found that NPDFs with support over the 
whole real axis, like the normal distribution, may introduce 
unphysical behavior and wrongly imply that values of <p < 0 
are accessible to A f (ip). For the lognormal and gamma dis- 
tributions it was found that <p — 1 defines the characteristic 
ice nucleation properties of the aerosol population. Although 
this was shown for deposition ice nucleation, it is likely that 
immersion and condensation freezing display the same be- 
havior as Eq. (36) is common to different heterogeneous nu- 
cleation modes. 

Singular and stochastic behavior were reproduced using 
the new formulation of N C (S U T,fi i... n ). In deposition ice 
nucleation S \ , 9, u, 7, ,v p , and a v play a role defining /f, 
however highly efficient IN ( 9 < 10°) are only sensitive to 
variation in 5), and 9. This is however not the case for less 
efficient IN for which changes in the cooling rate (expressed 
through the dependency of /f on u ) and variability in surface 
characteristics (expressed through the dependency of /f on 
cfy) impact /f. Thus, ice nucleation on IN with high 0 show 
features commonly identified with stochastic behavior while 
the ice nucleation spectrum of efficient IN (low 9) tends be 
consistent with singular behavior. This implies that ice nu- 
cleation on materials with multimodal NPDF would likely 
exhibit both, singular and stochastic characteristics. 

It was found that even if N c (Si, T,pi\ n ) is affected by 
variation in the cooling rate it may still be insensitive to 
temporal effects. This is because temporal effects may be 
masked by surface variability and by high or low nucleation 
rates at the experimental conditions employed. However 
high surface variability will not always lead to weak temporal 
effects as they also depend on the 5) trajectory followed in an 
experiment. It was also shown that even if temporal effects 
on N c (Si , 7, pu i ... n ) are significant, N c (S t , 7, pn... D ) may still 
be insensitive to variation in the cooling rate. Both cooling 
rate independence and constancy of /f over time were repro- 
duced using a time-dependent, stochastic approach. Thus, 
this behavior cannot be considered prove of the singular hy- 
pothesis and measurement of /f alone cannot discriminate 
between singular and stochastic ice nucleation. 

The temporal dependency of f{ shown in Fig. 4 suggested 
that IN measurements in cloud chambers with short residence 
time (~ 10 — 20 s) would tend to underestimate the IN con- 
centration. However it must be noticed that in this work it 
was assumed that active sites can be represented as adsorp- 
tion sites with high Thet- Nucleation on active sites must 
admit a time-dependent representation, however there is no 
guarantee that it can be approximated through Then Thus 
the temporal dependency of /f may differ from the results of 
Fig. 4 if active sites that cannot be represented as adsorption 
surfaces are present in the IN population. 

Comparison of derived N c (Si, 7, // | n ) against empirical 
data (Mohler et al., 2006; Welti et ah, 2009) showed that the 
NPDF of simple dust species like Kaolinite and Montmo- 


rillonite can be accurately represented by bimodal NPDFs 
whereas complex mixtures like Arizona Test Dust and Sa- 
haran Dust are well represented by three-modal NPDFs. It 
is however remarkable that for all of the dust species stud- 
ied these modes are consistently located around 9 1 ~ 10 °, 
§2 ^ 15°, and 9 3 ~ 20°, suggesting that the NPDF of dust 
from different sources can be represented by a linear com- 
bination of few lognormal modes. This feature was used 
to generate a parameterization of deposition ice nucleation 
on “generic dust”. Further measurements of N c (Si , 7, yu. 1 ... n ) 
covering the full range of /f (i.e. /f = [0 — 1 ]) are required to 
elucidate the fundamental reasons behind such behavior. 

The nucleation spectrum of soot was also investigated. It 
was found that deposition ice nucleation on low organic con- 
tent soot is well represented using a single-mode lognormal 
NPDF with 9 — 34.2°. The limiting /f of 25 % found by 
Crawford et al. (2011) was in good agreement with the pre- 
dictions of Eq. (53). A parameterization was proposed for 
usage in atmospheric models, however further experimental 
measurements are required to better represent the ice nucle- 
ation properties of soot. 

The theory presented here suggests that inferring the 
aerosol ice nucleation properties from measurement of 
5), onset at f{ ~ 0.01 may carry significant error. First, Fig. 3 
shows that 5), onset is strongly impacted by a v so that it may 
vary between samples of the same material or between freez- 
ing cycles. Second, .S) onse t represents the ice nucleation 
properties associated with the most efficient mode in the 
NPDF, which may be the less abundant (Table 2). Finally 
5), onset does not provide information on a ^ which may im- 
pact the inferred 9 (Fig. 5). These effects can be reduced by 
measuring Aj me d at f{ ~ 0.5 (corresponding to ip ss 1). Ide- 
ally, it is best to determine the /f vs. 5; curve which would 
also indicate the existence of multiple modes in the NPDF. 
For this, it is also important to measure the aerosol number 
area distribution instead of the number size distribution as it 
is the former, not the latter, what determines <p and hence 
N c (Si, 7,|t[ n ). 

Although this work focuses on heterogeneous nucleation 
in the deposition mode, it is readily suited for other nucle- 
ation modes. Heterogeneous freezing in the immersion and 
condensation modes share many features with deposition ice 
nucleation. Still, water adsorption on the particle surface 
(e.g. Seisel et ah, 2005; Kumar et ah, 2009), the effect of 
the solute on the interfacial tension, and the droplet size dis- 
tribution may play a role in defining fi. This will be detailed 
in a companion study. 

In a general sense, the NPDF can be understood as the 
combined distribution of ice nucleation efficiency and the 
area over which that efficiency maintains. Although the 
NPDF is fundamentally different from the contact angle and 
the active site distributions, they can be recovered from the 
NPDF when either the contact angle or the active site area is 
constant in the aerosol population. Thus using the theory pre- 
sented here the results of previous mechanistic approaches 
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can be understood within a more comprehensive framework, 
which at the same time is casted in a simple mathemati- 
cal form requiring fewer parameters, and unlike fits to ex- 
perimental results assigns precise physical meaning to such 
parameters. This work provides a physically-based method 
to reconcile theory, laboratory measurements and field cam- 
paign data, within a unified ice nucleation framework suit- 
able for atmospheric modeling studies. 
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